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A recent article suggested that the saturation of low energy spectral weight observed by X-ray 
absorption spectroscopy in the cuprates at high hole doping could be explained within the single- 
band Hubbard model. We show that this result is an artifact of inappropriate integration limits. 



Recent measurements [H 0] by X-ray absorption spec- 
troscopy (XAS) have shown an abrupt saturation in the 
oxygen K edge low energy spectral weight at a hole 
doping around x = 0.22 in several different overdoped 
cuprate systems. It has been argued [2j that this could 
not be explained within the single-band Hubbard ap- 
proaches commonly used to model the cuprates, possibly 
associated with the inapplicability in that doping regime 
of the basis state underpinning such approaches. Cal- 
culations have since indicated that the conclusions may 
be more general and more significant being appar- 
ently at odds with the three-band Hubbard model which 
the single-band model is intended to approximate. How- 
ever, calculations have recently been presented contend- 
ing that the saturation can be explained within a single- 
band model We examine this assertion. 

Since the U = 6t = 0.75W (where t is the nearest 
neighbor hopping integral and W is the band width in 
a 2 dimensional square lattice) spectra that Phillips and 
Jarrell present are, as they mention, continuous, with no 
gap below the upper Hubbard band (UHB), a nontrivial 
choice of upper integration limit is required to generate 
the integrated density of states shown in the inset to their 
Figure 1. The cutoffs used, in units of At, were lo c = 0.25 
for U = 6t = 0.75W at all dopings and oj c = 0.17 for 
U = 8t = W at all dopings, or oj c ~ 2J (J being the an- 
tiferromagnetic nearest neighbor exchange interaction). 
These somewhat arbitrary choices of integration cutoff 
are problematic to the point of effecting the qualitative 
doping dependence presented in their inset, for two rea- 
sons. 

First, these narrow integration windows fail to capture 
a large fraction of the low energy spectral weight (i.e. 
the weight below the upper Hubbard band) . Indeed, this 
cutoff is near the peak of one of the features of interest 
at low dopings, inherently missing around half of that 
peak's weight. Clear evidence of the appropriateness of 
the authors' cutoff can be seen in the unphysically low 
integrated weights in their inset — before even consider- 



ing dynamical spectal weight transfer (DSWT), which is 
central to their argument, slopes at low doping in a Hub- 
bard model must exceed unity per spin state to be phys- 
ically meaningful. DSWT should increase this slope sub- 
stantially at low dopings, particularly for the parameters 
chosen, which will tend to strongly enhance the spectral 
weight transfer. The slope in the inset for U — 8t = W is 
clearly well below unity, making those points unphysical 
and in disagreement with the results in their references 
[7] [3 and [8] [|. The slope for U = 6t = 0.75W is also 
below unity for all segments present (this curve extrap- 
olates to a positive intercept, not zero as in the physical 
system which is an insulator, because the bandwidth ex- 
ceeds U). 

A paper by Meinders et al [7j provides a useful demon- 
stration of the behavior expected within the Hubbard 
model. Meinders' low energy spectral weight increases 
strongly with doping for low dopings, as seen in their 
Figure 3. Meinders' points must be reduced by a fac- 
tor of two for comparison to those in [4[, which are per 
spin. This results in expected slopes of about 1.5 — 2 
for the most comparable parameters, well above unity. 
Phillips and Jarrell used a more sophisticated approach 
[HIH and in two dimensions rather than one, but the gen- 
eral arguments for spectral weight transfer are dimension- 
independent. The factor of ~2 disagreement at low dop- 
ings between Meinders and Phillips is anomalously large, 
and that the parameters chosen by Phillips ought to re- 
sult in even steeper slopes at low doping makes this dis- 
agreement particularly striking. 

As for the choice of an w c , apart from the reasoning 
that integrating the low energy spectral weight (i.e. the 
weight below the UHB) requires integrating all of the 
weight below the UHB, comparisons to existing theory 
0, 0] or experiment Q require integrating to the gap or 
minimum below the upper Hubbard band as was done 
in that previous work. Choosing integration cutoffs con- 
sistent with experiment is particularly important in this 
case, since the purpose of Phillips' paper is to explain 



2 



the experimental data. The integration cutoff must be 
in the broad minimum found well above [4| 's integration 
cutoff. The minimum below the UHB should be signif- 
icantly clearer in the U = 8t = W simulations, which 
were not shown. 

Second, and more crucially, the cutoff must change 
with doping. Since hole doping (removing electrons) 
manifestly shifts the chemical potential, to which all ener- 
gies in [ij are referenced, the low energy spectral weight 
and UHB are shifting with respect to the energy axis 
and integration window. This effect is clearly visible in 
Phillips' paper itself, where the UHB weakens and moves 
to significantly higher energy on doping, consistent with 
the results in Meinders' Figure 2 [7(. Indeed, Phillips' 
UHB shifts by roughly three times the width of the in- 
tegration window between x — 0.05 and x = 0.20; the 
peak of the low energy spectral weight shifts by about 
half that, but it should be gaining states on its low en- 
ergy side, so this is unsurprising. 

Experimental data [2J clearly show the addition of 
weight on the low energy side of the lowest-energy pre- 
peak as holes are added, moving the chemical poten- 
tial further below the UHB. The upper integration limit 
in this paper, located in the minimum between the low 
energy prepeak and the UHB as mentioned above, was 
not changed with doping because zero in this case corre- 
sponds not to the chemical potential but to the energy of 
removing an electron from a Is orbital, and should have 
negligible doping dependence (this electron does not ex- 
ist in the single-band picture). It should be noted that 
the experimental data are broadened by about 0.5eV due 
to resolution and core hole lifetime effects, masking small 
features such as any pseudogap; the fact that the UHB 
can be clearly distinguished experimentally from the low- 
energy weight, especially at low doping, is indicative of a 
correlation-based gap significantly larger than that used 
by Philips. Meinders Q had the benefit of a gap between 
the low energy spectral weight and the UHB, and inte- 
grated the former in its entirety and none of the latter. 
In this paper's second figure, a shift in the density of 
states relative to the chemical potential may be clearly 
observed. 

Failure to change the integration cutoff with doping 
can create Phillips' reported doping dependence as an ar- 
tifact, since the integration window effectively moves to 
exclude more spectral weight on doping. We reanalyzed 
Phillips' and Jarrells' spectra, integrating from zero to 
the minimum immediately below the UHB, using sev- 
eral different schemes to choose a cutoff (Figure [IJ . This 
captured roughly twice as much weight, gave physically 
reasonable slopes at low doping, and consistently failed 
to produce saturation at high dopings. The qualitative 
results were not sensitive to reasonable choices of cutoffs 
here, so long as the cutoffs were consistent. The reanal- 
ysed results have essentially the same doping evolution 
as that obtained by ID exact diagonalization [7[, and are 
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FIG. 1: a)Spectra from [J, reproduced for reference. 
b)Integrated low energy spectral weight per spin. Points from 
Phillips' inset may be compared with the results obtained 
when an upper integration limit near the minimum between 
the UHB and low energy spectral weight is chosen by five dif- 
ferent schemes. Curves reliant upon the position of the UHB 
are slightly suppressed at the highest doping due to the UHB 
becoming weak and indistinct. The most directly compara- 
ble curve from Meinders' Figure 3 [?J has been included — 
Phillips' calculations are for a very low U/W (W the band- 
width), and a choice of similar parameters by Meinders would 
have produced a curve above that shown. The bandwidth in 
ID (Meinders) is W = At and in 2D (Phillips) is W = 8t. A 
line of slope 1, the minimum physically reasonable slope at 
low doping in a Hubbard model, has been included for refer- 
ence. The shape of the curve is not sensitive to the choice of 
cutoff near the minimum. 

just as inconsistent with the experimental data. 

A further qualitative explanation for the saturation is 
offered in Q's final paragraph, in which the closing of 
the pseudogap halts the dynamical spectal weight trans- 
fer (DSWT). In fact, this is seen in the reanalyzed weights 
(Figure [1]), where the slope returns to roughly unity 
around x = 0.2. Phillips' and Jarrell's weights do not 
return to a slope of unity when the slope ceases to be 
enhanced by DSWT, they fall to a slope of nearly zero. 
The slope per spin in a Hubbard model without DSWT 



3 



is unity and with DSWT is well above unity, at least at 
low dopings, and in a Fermi liquid is 0.5. The Fermi 
liquid slope is easy to see: If we start electron doping a 
completely empty (x = 1) Hubbard system, then each 
electron blocks one electron addition state; the UHB has 
vanishing weight here, so the slope is 1 per electron, or 
0.5 per spin. If we had a U — system in this picture, we 
would start from an empty (x = 1, not half filled) band 
and add electrons at a slope of -1 (-0.5 per spin). Since 
there are two electron states per Cu (or one per spin), 
our initial value would be two (one) and we would cross 
half filling at a weight of 1 (0.5). In a Hubbard system 
at half filling there is no low-energy weight, and remov- 
ing one electron creates two low energy electron-addition 
states, the second of which had previously been blocked 
by the electron thanks to the on-site repulsion U. This 
leads to a slope of 2, or 1 per spin. Exact diagonalization 
in ID has shown a smooth evolution from one regime to 
the other, with a slope of 0.5 per spin at x = 1. Having 
a slope of zero would require that one of these regimes 
evolve into the other nontrivially via a saddle compatible 
with neither picture. 

Having reanalyzed Phillips' and Jarrell's calculations 
using integration limits consistent with those used on the 
experimental data and in previous theoretical work, we 
have shown their results to be consistent with ID exact 
diagonalization and inconsistent with the experimental 
data. Unfortunately, we are forced to conclude that the 
substance of [3| is based entirely on an artifact of an 
inappropriate choice of integration cutoff. There remains 
no evidence that single-band Hubbard approaches can 
explain the XAS data at high dopings. 
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Note: in the original PRL Q, the upper integration 
cutoff for Tl-2201 was incorrectly reported as 529. 4eV, 
the same as for LSCO. In fact, the upper integration 
cutoff used for Tl-2201 was 528.8eV. All other cutoffs are 
correct in the published paper. This would not have come 
to our attention were it not for the work of P. Phillips 
and M. Jarrell, to whom the authors are grateful. 
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